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INTRODUCTION 

In  an  academic  sense  alone,  a  mathematical  investigation 
of  the  stability  of  towed  bodies  and  systems  is  of  interest 
only  if  it  adds  to  the  existing  knowledge  in  the  general  field 
of  the  course-keeping  qualities  of  vessels.   In  a  more  practi- 
cal sense,  it  is  an  important  and,  in  fact,  necessary  primary 
design  aid  which  yields  specific  criteria  to  be  satisfied  for 
stable  operation.,  These  criteria  are  composed  of  terms  which, 
for  the  most  part,  have  distinct  physical  significance  and 
which  can  be  varied  in  design  studies.  Through  an  appropriate 
analysis  of  such  terms,  desirable  operating  characteristics  may 
be  found  to  a  considerable  degree  of  accuracy  while  such  a  de- 
sign is  still  in  the  "drafting  board"  stage. 

This  paper  deals  with  the  analysis  of  a  system  composed 
of  a  towboat  or  tug  and  a  single  barge.   Both  bodies  are  con- 
sidered to  be  unrestrained,  hence  the  stability  of  the  system 
will  be  dependent  on  the  characteristics  of  the  component 
partso  A  simple  automatic  control  is  specified  for  the  towing 
body  in  order  that  the  system  should  have  directional  rather 
than  dynamic  stability  (or  instability). 

In  general  the  methods  followed  are  those  of  Euler,  Routh 
and  Hurwitz  and  deal  with  systems  of  simple  linear  differential 
equations  only„  Laplace  Transformations  were  not  applied, 
since  the  author  has  found  that  many  engineers  are  unfamiliar 
with  th-'  properties  of  this  concept. 
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The  subsequent  analysis  presupposes  the  follcfwing  basic 
simplifying  assumptions: 

1.  Both  bodies  have  freedom  in  the  horizontal  plane 
alone,  i.e.,  roll  and  pitch  are  ignored. 

2.  The  tow line  is  assumed  horizontal  and  taut  at  all 
times. 

3.  Angular  displacements  are  considered  to  be  small,  so 
that  second  order  and  higher  terms  may  be  neglected  in  the 
final  equations  of  motion. 

With  these  premises  in  mind,  consider  a  set  of  right 
hand  orthogonal  axes  fixed  relative  to  the  earth  and  denoted 
^y  i^  0  "^  o) "     Consider  also  similar  systems  of  Bxes-  {^,y)   and 
(x,y)  fixed  in  the  barge  and  towboat  respectively.   Such  a 
system  is  shown  in  Fig.  1  together  with  other  symbols  neces- 
sary to  define  completely  the  orientation  of  the  system  at  any 
time.  A  complete  and  detailed  nomenclature  is  given  in  Appen- 
dix T. 


Figure  1 
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Kinematic  Relations 

CertaJn  expressions  are  obtained  from  the  towing  geometry. 

These  are  presented  without  comment  based  on  Fig.  1  of  Ref.  1. 

Vo,''  J  V  Sin  fi//  -  ^)dt  >  x^  Sii»A  =  y^^  *    ?  Sin  (v  -  4') 

x^^=  j  VCos  (^  -  (i)dt  ^  x^  GO90  ^'op-   ^Cos  (r,-4>) 

y^^  =  J  F  Sin  (4>  -  P)dt  ^  x^  Sin0  (1) 

x^     =  J  P  Cos  C0  -  J^)dt  +  J    CO80 

.:  7}  =7<  +  arc  Cos  f"?  Cos   (yp  -  ^)dt      J  J  Cos   (^-/3)dt  +y''  Co80  -y  Cdb^j 

=  17'  +  arc  Sinrf|  Sin  (4>  -  0)dt  -  f  |  Sin  (J^-^)dt  +y  Sinu';     y'  Sin0 

^  -^=   rj-^P  (2) 

r^=TCos^         .       'ry=-TSinf      ,       N^^f-xJSin^ 

(3) 
Tj  =  -  r  Cost)     ,       Ty  =  T  SiuV  ,       /V^^=  z^  T  Sin^ 

Differential  Equations  of  Motion 

Since  axes  fixed  in  the  bodies  are  essential  to  the  de- 
velopment it  is  necessary  to  follow  Euler,  Ref.  2  in  determining 
appropriate  mass-acceleration  terms.   For  the  pui^jose  of  con- 
tinuity these  expressions  will  now  be  derived. 

Choose  a  set  of  axes  fixed  in  space  and  a  set  of  axes 
fixed  in  the  body  as  shown  in  Fig.  1.   If  P^   and  Py   are  the  mo- 
mentum expressions  relative  to  the  axes  fixed  in  the  body  then 
(according  to  Newton)  the  force  components  along  the  axes  fixed 
in  space  are  equal  to  the  components  of  the  time  rate  of  change 
of  momentum  along  these  axes»  Thus 


X    '  ^  p     =  ^ 

0    dt   "'    dt 

P,  QxW'- 

/'v  Sins- 

Y     ^    d    u       =    d 

^^   dt  ^'y   dt 

P,'  SirW/  * 

/'y  COS0 

also 

Xg  =  X  Cos^J  -  Y  Sin/' 

Y^  -  X  SirW^  + 

y  rm0 

•••  X=^-Py^ 

Y  =Py-  P^4' 

where 

Pj,    =  M,    U           .            Py    =  M,    V 

• 

Then 

X  =  V/j  u    -  Wj  V  «// 

y  =  Wj  0  ^  Wj  u  v/) 

Similarly 

^  =  Wj  U  -  My  v^ 

y  =  Wj  V  ^  Wj  n  0 


^  =  ^B  ^  ^  '^ 


(6) 


(7) 


(6) 


(9) 


The  moment  of  inertia  terms  are 

N  =  I  ^        and   /V  =  7  0  (10) 

We  consider  here  Ref ,  1  and  Ref.  2  and  note  that 
Wj  =  W  Ci  +  kj)        .        M,  =  M  (1   +  fej) 


(11) 


where  M   is  the  absolute  mass  of  the  ship;  Jg   is  the  moment  of 
inertia  of  the  body  about  a  vertical  axis  taken  through  its 
center  of  gravity,  I^   is  the  moment  of  inertia  of  the  fluid 
displaced  by  the  body  taken  about  a  vertical  axis  through  its 
center  of  gravity  and  k ^,    k^   and  k     are  appropriate  coeffi- 


cients  of  accession  to  inertia.   Similar  expressions  may  be 

written  for  the  "bar"  terms,  thus 

r  =  W  Ci  ^  fe,)      M.^  =  M  (1  *  k  ) 

(12) 

We  are  now  in  a  position  to  write  the  differential  equa- 
tions of  motion: 

1.  Barge 

a)  Longitudinal  Force 

Mj  u  -  Wj  u  J/  =  A'  (u,v,i)   +  T^  (13) 

b)  Lateral  Force 

Mj  i  +  Mj  u  di  =  Y  (u.v,i>)   +  Ty  (lU) 

c)  Yawing  Moment 

JV^  =/V  (u.v^)  +  Nfj^  (15) 

2.  Towboat 

a)  Longitudinal  Force 

Jijii-M^vii  =  X(u.v,^,'l)*f^  (16) 

b)  Lateral  Force 

W,  v  +  W,  ii  0  =  F  fu.v.^J;  +  fy  (17) 

c)  Yawing  Moment 

70  =/V  (xx.v.Jj^)   +  /V^T-,  (18) 

3.  The  towboat  rudder  equations  for  rudder  angle  pro- 
portional to  heading  change  are: 

a)  Rudder  order 

^*  =  a„;/;  (19) 

b)  Constant  time  lag 

dt 
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Linear  Simplifications 

Unles:-  the  previous  systems  are  linearized,  a  classical 
solution  iu  •■xii  umely  difficult  to  obtain,  in  fact  in  a  non- 
linear form  the  problem  can  best  be  treated  through  the  use  of 
automatic  computing  devices.   However,  only  a  closed  form  solu- 
tion will  be  considered  herein,  hence  the  systems  will  now  be 
linearized. 

It  is  first  assumed  that  K  ^  7  ,  i.e.  that  both  the  tow- 
boat  and  barge  are  traveling  at  the  same  speed  and  that  both  ^ 
and  I^  take  on  only  small  values.  Then  the  following  relations 
are  obtained  from  Figure  1 

u  =  K  Cos/^  :z     y 

V  =  -y   Sin^  :^  -y[^ 

a  =     y   Cos/3  ;i  1/  •  ^      ' 

V  =  -K  Sin/T~  -V)F 

Noting  these  relations,  the  truth  of  the  following  iden- 
tities is  obvious 

A  (n.v.^'h)         =  X    (l\C') 

Y  (u.v,yh)        =   Y    (b,^Jj)  (22) 
N    (u,v,u<)         =  /V  (P.'^!l) 

X  (u.v.i),!))  =  X  (B.^.B) 

Y  (u,v,i,h)   =  Y   (B,$.J^) 

/V  (u,v,d!,J)    =  N   C/^>,») 

Limiting  all  displacements  to  the  first  order  of  small 
quantities,  relations  (22)  may  be  expanded  in  Taylor's  Series 
to  obtain  the  following: 
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fid^.H')  ^  ^  '  1^  <'^  -  ^e^  ^  U  <''^'  -  '^e^  (23) 

3/3  bA"  Bd 

—   -T  -  -  _    3/V-—     '^N     —      —         3^__ 

where  the  subscript  e  denotes  the  equilibrium  values  about  which 
the  expansion  takes  place.   Since  it  is  assumed  that  initially 
the  system  is  on  straight  course  and  in  line  it  is  obvious  that: 

and 

/3,  =  0e  =  ^e  =  ^e  =  ^  =  ^   • 

In  addition,  for  reasons  of  symmetry,  Ref.  3, 

3/3   BB"   30   cKi7   BT 
Relations  (22)  then  reduce  to: 


/V   =  /V^/3  +  /V^>A 


X   =  X^  Y  -  Ys^  ^   Yph  ^   K-jS 


N   =  N^0  ^  N^4'  +  /VjS 


(21.) 


where  the  subscript  notation  has  replaced  the  conventional  par- 
tial derivative  form. 
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Continuing  the  linearization  to  include  expressions  (1), 
(2),  and  (3)  the  following  relations  are  established 


L  -T 


T,--i 


(25) 


p  ,       p  - 


^    (u'-^)Jt    +  y  </'  -y  « 


7y  ^''l'"   {j  ('^-^)^'   -Jj 


^n  =  ^/ 


(26) 


In  a  similar  fashion  expressions  (O),  (9)  and  (10)  become 

V  =  0 

^f  =  /  iv; 

K  =  0 

Y  =  ¥^  V  J'  -M^V  0 

/V  =  J^ 

With  these  linearizations  the  equations  of  motion  become: 

1.   Barge 

a)  Longitudinal  Force 
0  =  Vg  +  T 


b)  Lateral  Force 


(27) 

=  Y/  +  Y^4> 

(28) 

-T 

t'.  +  r^(w-e)dt  - 

Cj-    (v'-^)dt 

-T 

~  X  „                   X-    — 

Y .  -  jp ._ 
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c)  Yawing  Moment 

-x^,r[fp.y -^P^]  (29) 

2.  Towboat 

a)  Longitudinal  Force 

f^-\-T  (30) 

b)  Lateral  Force 

-tI^  (4'-^)dt  *fp:p-jP^]    (31) 

c)  Yawing  Moment 

1$  =  N^^N^^^NiJ  -  ipTW^^Jj(yP-0)dt  -  jj  (4^-Ii)dt'] 

^^pT^jP^'-jP^j  (32) 

d)  Rudder  Order 

^'  =%^  (33) 

e)  Constant  Time  Lag 

r^i.tf  (3U) 

at 
Equations  (27)  and  (JO)  represent  steady-state  conditions 

and  need  not  be  considered  further  in  the  stability  analysis. 

Non-Dimensionalization  of  Equations  of  Motion 

Since  the  hydrodynamic  forces  may  best  be  expressed  in 
terms  of  dimensionless  coefficients,  one  divides  equation 

(26)byffV         5   (29)byf/V2 

(31)  by|/V  :    and   (32)   by  |/V2 

Furthermore  the  following  transformation  is  introduced  for  the 

Vt 

independent  variable  t  ,     s   — — 

Here  s  is  the  number  of  tow  lengths  traveled,  t   the  true 

time,  L   the  tow  length  and  V  the  linear  velocity  of  the  system. 

Since  V  is  considered  substantially  constant  the  following 

relation  is  written  v  "^  ,  •  (35) 

at      L 


thus  the  following  expressions  are  obtained: 


_  ir,  _ 


ds 

ds 

dt 

» 

0 

V 

0    =   # 
ds 

(is 

(if 

V 
L 

7 

i 

^    dfi      ds     ^    f/ 
ds       dt 

hi 

ds 

ds 

dt 

= 

•  1 

0 

V" 

0 

ds        dt 

> 
L 

2 

t 

-  dB     ds    ^   fl 
ds      dt 

In  addition 

■i      =   . 

r 

2 

-r 

' 

n 

'       _          2 

2^ 

' 

n      = 

I 

r-i    I'-i     I'-i 


(36) 


_  / 

2>t 

Similarly: 

y'. 

2^ 

y'i 

2^^ 


2-^ 


N 


N 


0 


2^ 


T     = 


iry' 


/v;  = 


/v. 


2 


And: 


9  /*•  2 


)9 


Y  -  =        N' 


Pi    2y2 


ir^^ 


H 


Using  the  foregoing  relations  the  equations  of  mot..   .- 
come: 

1 .  Barge 

a)  Lateral  Force 

<i'-  --: /^'  -->'/  I'V  (37) 

b)  Yawing  Moment 

n/'^'  -/v'/  -/v'^j;^  (36) 

2 .  Towboat 

a)  Lateral  Force 

<  l' i'  -^2  l' ^' -y'-^-y'i^'  -y'-^         (39) 

b)  Yawing  Moment 

—  /    '  r- 

c)  Rudder  Order 

8'  =  a^  ^  (Ul) 

d)  Constant  Time  Lag 

I'  =  6  .  5  ^  (Ii2) 

as 
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STABILITY  ANALYSIS 

The  equations  of  motion  of  the  bodies  under  consideration 
together  with  the  coupling  effects  due  to  the  towline  con- 
straint and  the  imposition  of  the  control  functions  on  the  tow- 
ing body  have  been  presented  in  the  previous  section. 

These  equations  are  of  general  significance  within  the 
limits  of  linear  theory  with  the  exception  that  a  constant  time 
lag  has  been  assumed  for  the  control  function.  This  assumption 
has  been  made  purely  to  show  that  time  lags  of  a  sort  can  be 
included  in  the  development.   However  in  the  subsequent  analy- 
sis the  control  time  lag  «  will  be  assumed  to  be  zero.   No 
loss  in  generality  will  occur  since  it  is  evident  that  the  in- 
clusion of  such  a  term  will  only  serve  to  increase  the  order 
of  the  stability  equation  by  one.   Its  exclusion  will  not  make 
the  analysis  any  less  valid  but  will  to  some  degree  serve  to 
simplify  it. 

To  remove  the  integral  signs  in  equations  (37),  (38), 
(39)  and  (ijO),  we  operate  on  them  by  -r     and  hence  raise 
their  order  by  one.  The  equations  then  become: 

1.   Barge 

a)   Lateral  Force 

.,7'-  y'iY^fY'',  -;']-' 'f'^ 


/  A' 


,  /  -/ 


-W ^'  -y'/-f^-fl^p  ^ 


-  1   0+X  e  =  0 


(U3) 


b)  Yawing  Moment 


r         ^      r  1' 

r  r 


^'    j'     ^1 


4> 


2. 

Towboat 

a)  Lateral  Force 

\l     -' 

-./iVt'-7Vl'-i>-i7',;  J' 

t'           t'         - 

.  1 

1  — 

r     r      ' 

b)   Yawing  Moment 


(Uh) 


(U5) 


?' 


?' 


^' 


—    4^ 


(U6) 


x'    T'             x'    t'  -I-' 

-     P         ^>     +  _£ /'-/Vt5  =  0 


?' 


?' 


V 


c )  Rudder 
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5  =  a„  li/ 
o 


(U7) 


These  are  linear  homogeneous  equations  in  the  dependent 
variables  ^  ,  /^  ,  0  ,  ^  ,  ^ ■    hence  their  solutions,  excluding 
the  highly  improbable  cases  where  the  roots  are  equal,  are 
given  by  the  following  equations: 


v//  =  I  0^  « 


oris 


I 


<ri5 


i 


(US) 


6=1     6.  e'^^^ 


.U4- 


It  is  our  aim  to  determine: 

(1)  the  condition  for  which  nontrivial  solutions  for  the 
value  7j  exist  in  the  system 

(2)  the  condition  that  these  <^i   values  are  such  that  the 
motion  be  stable,  i.e.,  that  the  real  parts  of  all  the  <^i  values 
be  less  than  zero. 

In  order  to  fulfill  condition  (1)  above,  we  substitute 
expressions  (U8)  into  equations  (lt3)  through  (U7)  and  realizing 
that  the  principle  of  superposition  of  solutions  holds  for 
linear  differential  equations,  drop  the  sujrnnation  notation  and 
obtain  the  following  relations: 

1.   Barge 

a)  Lateral  Force 

([•■7-  y'iy^^jY'^  *?']^.^f.]^. 

b)  Yawing  Moment 


( 


r 


I'       '     r 


^"12   /\  =  0  .  (50) 

r 


2.   Towboat 

a)  Lateral  Force 
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b)  Yawing  Moment 


(51) 


—  It'  I 


n-'/^-.^-  ^;^^,^*:^r/x-;_?';.,  .ilb 


-('vW. 


.!^n-.r-;  y/^;  .^^ 


r 


r 


r 


'^-^ ^-f^\-.h-  ^ 


c)  Rudder 


b     =  a     ^Jj  ■ 
I  o  ■  I 


(52) 


(53) 


We  now  have  a  system  of  linear  homogeneous  equations  in 
^i*  ^i  f  ^i  *  ^'i  >   ^^'^   ^i'     '^^^   condition  that  a  nontrivial  solu- 
tion exists  among  them  is  that  the  determinant  of  their  coef- 
ficients be  identically  equal  to  zero.  This  means  that 


^11 


^21 


*31 
^51 


'12 


22 


*32 
'S2 


^13 


^23 
'33 

'a3 

'53 


'3U 

u 

'5U 


'15 
'25 
'35 
\S 
'55 


=  0 


(5U) 


where: 


(55) 
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Ktr'        -L 


^13 


^Ih 


a .     +     ;/ 


— (7        *      ■) 

r     '    ^ 


I 

r 


(55) 


^15 

^23 
"2U 
^25 

^32 

"33 
"3U 
"35 

"Ul 
"U2 


x'      T' 


x'    r' 


'^jllz'.,  .!iL_''' 


;' 


r 


=     0 


r     I 


/' t'        t'' 


-r-^.^f 


1 

2' 


-y  j^^ 


jc„      X  ^  Jl.  T  X ^      T 


r 


-.  --^ 


r 


r 


(56) 


(57) 


(58) 


%} 

=  ("'/'= -.'-n-.'^V^' 

'I'-'r-f 

■        r 

%k 

r-'       >„'  ''1 

%s 

=  -(V|   c. 
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^52  '  ^  (59) 

"5U  =  '         ^?^  =  -1 

Expan:;ion  of  this  determinant  leads  to  the  following  rela- 
tion: 

j=o  J  (60) 

where : 

Cq  =  (ah  -  id)(AH  -  LD) 

C^  =  (ah  -  >^d)(AJ  +  BH-LF-MD)  +  (aj  +  bh -/f  -  md)  (AH  -  LD) 

C2  =  (ah  -  /d)(BJ  +  EH  ♦  LG  -  MF  -  ND) 

+  (aj  +  bh  -^f  -  md)(AJ  +  BH  -  LF  -  MD) 

+  (bj  +  eh  -  ^g  -  mf  -  jg)(AH  -  LD) 

C^  =  (ah  -  i  d)(EJ  *  Mg  -  NF) 

+  (aj  +  bh  -  /f  -  md)(BJ  +  EH  +  Lg  -  MF  -  ND) 

+  (bj  +  eh  -Jg  -  mf  -  jg)  (AJ  +  BH  -  LF  -  MD)    (6I) 

+  (ej  -  mg  -  jf)(AH  -  LD) 

C^  -  (ah  -  i  d)  Ng  +  (aj  +  bh  -  /  f  -  md)(EJ  +  Mg  -  NF) 

+  (bj  +  eh  -  ^  g  -  mf  -  jg)(BJ  +  EH  +  Lg  -  MF-ND) 

+  (ej  -  mg  -  jf)(AJ  +  BH  -  LF  -  MD) 

-  jg  (AH  -  LD)  ^  g^£  L 

C^  =  (aj  ■»■  bh  -/  f  -  md)  Ng 

+  (bj  +  eh  -  ^  g  -  mf  -  jg)(EJ  +  Mg  -  NF) 

*   (ej  -  mg  -  jf)(BJ  +EH  +  Lg  -  MF  -  ND) 

-  jg  (AJ  +  BH  -  LF  -MD)  +  {g^ J   M  +  g^  mL) 
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C^  -  (bj  *  eh  -  ^g  -  mf  -  Jg)  Ng 

♦  (ej  -  mg  -  jf)(EJ  ♦  Mg  -  WF) 
-  jg  (BJ  +  EH  +  Lg  -  MF  -  ND) 

*   (g"i  N  +  g^mM  +  g^j  L)  (61) 

C^  =  (ej  -  mg  -  jf)  Ng  -  jg  (EJ  +  Mg  -  NF) 
+  (g^mN  +  g^j  M) 


'8 
where: 


Cq  =  (-jg^N  +  gjM) 


d  =  -'"2  /' 

f  =  -y^^  .       (62) 

g     ^. 

A  -  -'"2  ^p  / 

E  =  -[r'  *  r^^  a„  r'j] 
D  ^-s;^ 
F  =  -F'^- 


H  =  "2'  V/ 
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L-[n7'^  >  i^;  -^;'  I'']  (63) 

N  =  [r^  ^V-  '^V  *  '^o  ^-V  YJ-    Nj)] 
To  fulfill  condition  (2)  on  page  (iLi),  we  need  look  to 
the  theory  of  equations  as  developed  for  polynomials.   In  Ref, 
3,  we  find  that  the  conditions  for  stability  (the  real  parts 
of  all  the  roots  are  less  than  zero)  are  stipulated  by  the 
following  considerations: 


(1)  all  C.  >  0 

J 

(2)  all  determintal  relations  obtained  from 


0 
0 


0 
0 


0 
0 
0 


0 
0 
0 


0 

^8 

s 

^6 

c 

0 

0 

0 

^8 

c 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

By  striking  out: 

first  the  first  two  rows  and  columns 

>  0 


0 


I.e. 


0 


'3     "2 
second  the  first  three  rows  and  columns 


I.e. 


'0 


0 


0 


(6k) 


>   0 


(65) 


(66) 


(67) 
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and  so  on  until  all  such  expressions  have  been  vn'itten  in- 
cluding (6^). 

Tne;;e  inequaiiLies,  when  satisfied,  inaicate  si-ability 
for  the  system  under  consideration.  For  any  two  particular 
bodies,  a  towline  length  and  towpoint  locations  on  both  the 
barge  and  towboat  may  be  determined  to  give  desirable  stability 
characteristics.  This,  however,  requires  considerable  numeri- 
cal analysis  as  well  as  time,  unless  suitable  computing  fa- 
cilities are  available. 

The  development  of  the  theory  as  presented  here  is  con- 
sidered to  be  a  basis  for  future  numerical  investigation  of 
the  directional  stability  of  any  two  bodies  separated  by  a  tow- 
line  of  arbitrary  length  attached  to  arbitrary  points  along 
the  centerline  of  each  body. 
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A  SIMPLIFIED  THEORY 

The  foregoing  analysis  has  been  devoted  to  the  direction- 
al stability  of  a  towed  configuration  of  two  bodies,  separated 

by  a  towline  of  arbitrary  length,  attached  at  arbitrary  points 
along  the  center line  of  each  of  the  bodies,  and  further  con- 
fined by  a  control  rule  for  the  towing  body  proportional  to 
its  heading  deviation. 

It  is  evident  from  this  analysis  that  the  results  ob- 
tained defy,  without  considerable  numerical  calculation,   ra- 
tional qualitative  results.  Thus  in  order  to  correlate  quali- 
tative as  well  as  some  quantitative  predictions  with  data  from 
test  measurements,  certain  rather  broad  simplifications  will 
be  made. 

It  will  be  assumed  that 

(1)  the  towline  is  horizontal,  and  taut  (as  before); 

(2)  the  barge  does  not  roll  or  pitch 

(3)  the  towboat  proceeds  on  a  straight  and  prescribed 
course  regardless  of  the  motion  of  the  barge. 

With  this  latter  assumption,  all  the  towboat  functions 
vanish  and  only  the  barge  functions  remain. 

It  is  also  convenient  to  define  the  transformation  and 


dimensionless  coefficients  presented  on  pages  9  and  10  as: 

(68) 


.  =^ 
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•P  -  •!>' 


'      if 


p     eD'^v' 


,  I 


0  '  0 


V 

I 


N 


'2  p    #3 

No 


fl       Ep'^y'' 


^..'[i\ 


U 


I  T 

T     = 


£.0    2y7 


(66) 


Y  .;,  = 


<P        P-  6  '^y 


N 


'A 


"/»        P  0  *Y 


f/ 


With  these  in  mind,  equations  (ii3)  and  (Uii)  reduce  to: 
Lateral  Force 


■  ■  '        /    ■  /     T  / 


(69) 


Yawing  Moment 


"'0-/v'is^"  -^  -;  t'[i  -  "-^]^' 


(70) 


Reasoning  similar  to  that  presented  on  pages  13  through  17 
leads  to  the  stability  equation 


j  =  0  J 


(71) 


where 


r       -        ^        ^ 
Cq  -  n   mj 
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The  motion  will  be  stable  if  and  only  if  all  the  C 


J 


values  are  positive  and  the  following  relations  are  satisfied: 


and 


<=! 

s 

^' 

S 

^2 

=1 

^=0 

0 

0 

^  = 

0 

=2 

'=1 
'=3 

c 

0 

0 

0 

c 

=  C^  C^  -  Cq  C3  >  0 


(73) 


(71.) 


In  the  case  of  the  untowed  body,  it  is  noted  (Ref .  1) 


that 


/  /  _  _ 


N  4,y  0*  ('^i-y  ^)  N 13=  -'^  "2  ^1  ^2 


(75) 


where  cr^   and  0^2  are  the  dynamic  stability  indices  of  the  un- 
towed body.  Then  the  C.  values  may  be  written  as 

r      -       '      ^ 
Uq  -  n  mj 

Cj  ■=  [»p  r'.;  ri  ♦  i£-)  *  -.'J'*  „' .,'  ?,  5,]  (76) 


'^h' 


'p  ^e-"^ 


I' 

r 


Without  further  mathematical  manipulation,  the  following 
analysis  may  be  made:  by  definition,  all  the  C  .  values  must 
be  positive.  Then: 

1)  Cq  ^  0  always 
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2)  C^  >  0  for: 

a)  all  dynamically  stable  untowed  bodies,  since 

a,  and  cr^  <  0  ; 

b)  all  dynamically  neutrally  stable  untowed  bodies, 
since 

o-j  =  0  and  o'j  <  0  ; 

c)  all  dynamically  unstable  untowed  bodies  for 
which 

1^1  I  <  1^2  I. where  o^j  >  0  and  ^j  <  0  . 
As  a  consequences  of  this,  a  towed  body  cannot 
be  stable  on  course  if   I^J  -*  1^2 1   ^^<^ 
1^1  I  >  0. 

3)  C2  >  0  for: 

a)  all  dynamically  stable  untowed  bodies,  since 
the  towing  functions  are  positive  and  since 
both  aj  and  ct^  are  negative; 

b)  all  dynamically  neutrally  stable  untowed  bodies, 
since 

0-^   =  0   and  CTj  <  0  ; 

c)  all  dynamically  unstable  bodies  for  which 


\x)     C-  >  0  for  all  towed  bodies  (see  item  5  below), 

since  by  definition,  {^p     ^  p  ~  ^  p)    ^0,  and  usu- 
ally (J''^  -  "i'   -^  "i  )  >  0,  /V';  >  0  and  T^and 
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5)  C,  >  0  by  definition.   Since 


^f  *''^-'vV 


%-"■ 


-N'p     >     0 


this  condition  may  be  fulfilled  either  by  increasing 
X       (moving  the  towpoint  farther  forward),  or.  if 
this  is  in  a  practical  sense  impossible,  by  increas- 
ing the  drag  in  such  a  way  that  the  static  moment 
rate  will  remain  unchanged  or  is  decreased.  The 
condition  C.  >  0  defines  the  location  of  the  tow- 
point:   it  must  always  be  a  distance  x'p    forward  of 
the  center  of  gravity  of  the  ship,  which  is  greater 
than  the  distance  between  the  center  of  gravity  and 
the  center  of  pressure  of  the  static  force  Y^0. 
Assuming  that  the  C .  values  are  all  positive,  it  is  ne- 
cessary  to  examine  the  determinantal  relationships 


-n'n/  faj+5j)  5-j  aj  -T'l^^aj+a^j  g  +  fij  >0    (77) 


where 


Q"p   <   i'*f>*f, 
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and 


and  Q.   ft,  and  n'  m^  must  be  positive.  The  following  cases 
will  now  be  examined: 

i)  Untowed  body,  dynamically  stable.  Here, 

always  for  a'j  a^  >  0  and  C^i  +  ^j-^  *^  0*  Thus,  the 
sense  of  the  inequality  depends  on  the  magnitude  of 
T  /?;  hence,  to  achieve  stability,  it  is  desirable  to 
keep  this  term  small.  This  may  be  accomplished  by 
proper  selection  of  x'     ,   the  towpoint,  and  /^  ,  the 
towline  length. 

ii)  Untowed  body,  neutrally  stable.  The  first  term 

vanishes  since  ^j  0^2  ~  0;  thus  ft  must  be  made  small 
in  order  to  satisfy  the  inequality,  and  insure  sta- 
bility. To  accomplish  this,  the  same  reasoning  is 
applied  as  in  a). 

iii)  Untowed  body,  dynamically  unstable.   Here,  ^j  >  0 

and  as  long  as  CTj  ^  ^  and  |ctj |  <  l^j |  then 

-nm^      (^i   +  ^2^  "^l  ^2  ^  0 
Thus,  for  stability,  it  is  desirable  to  make  ft  as 

small  as  possible  and  to  increase  the  drag  in  such  a 

fashion  that  the  inequality  is  satisfied. 
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b)   A, 


^  "i  "2  "3  -  "i    %  ■  "O  "i  '   ° 


(78) 


=1  % 


>      0 


If  it  is  assumed  that  the  conditions  specified  for  Aj 
have  been  me+ ,  then  the  remaining  factor  to  consider  is  the 
gnitude  of  C,   (C. /C^),   C,  is  determined  by  the  values  of 


ma, 


CTj,  CTj  and  n'n^     ,    and  0  <<  C, /C  <<  1.   It  is  therefore  de- 
sirable to  make  C,  as  small  as  possible,  still  keeping  it 
positive,  and  thus  insure  the  stability  of  the  system. 

The  criteria  of  directional  stability  will  now  be  deter- 
mined for  the  case  where  the  towline  has  increased  in  length 
to  a  point  where  the  terms  containing  its  reciprocal  are 
negligible.  Mathematically,  this  is  achieved  by  letting  the 
towline  length  increase  without  limit.   In  this  case,  it  is 
found  that 


^py'0-N'^y]>o 


Since  n\  m'    ,  T^   and  Cx  "^  Y 'q   -N^q)       are  positive,  this  ex- 


pression  will  be  positive  if 


l^m/  (a,   >  ¥^)(x^  t'^  n'a^,  a,)  -  (x^  Y '^  -    \>  T  '  <j 


or 


(79) 


CV  ^'0    -'^  V  ^''<  -"2"  (^l    "^2^  (^p      '^'^   "'^1  ^2^- 


I 
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For  all  dynamically  stable  ships,  the  right-hand  side  of  the  in- 
equality is  positive}  thus  (x^     Y  p  -  N  p)T    must  be  kept  as 
small  as  possible  if  the  resulting  motion  is  to  be  airection- 
ally  stable.  For  dynamically  unstable  ships,  the  directional 
stability  depends  on  (x       T    +n^  a,  CTj)  >  o.   Since  in  this 
case  5-j  CTj  <0,  it  must  be  specified  that  for  directional  sta- 
bility,  X   7  >  -  n  cTj  cTj  in  addition  to  the  condition  that 


(^p  ^'b   -^V  t'  <   -'^2     (^,  *^2)   (^p     l'-  na^  ^^) 


P 
p   ^   fi   ~"   /3'    '       "        '"2   '^1  ■  ^2'  "'p 

To  investigate  the  consequences  of  decreasing  the  towline 
length  to  zero,  i.e.,  either  towing  the  ship  from  its  bow  or 
with  a  bridle,  we  discard  the  terms  containing  the  towline 
length  ]     in  equations  (69)  and  (70)  on  page  22  since  they  have 
no  significance  in  this  situation.  The  resulting  stability 
equation  is  then  a  cubic  of  the  form  of 

(80) 


3-j 

1 

c. 
J 

a,  -0 

where 

So 

= 

n      m„ 

^1 

= 

("'y'^  - 

"2 



,..'           1 

/ 

fli^  N   ^)    ~    -n       m2   (^i    "*■  0-2) 

C3  -  rx;  Y'p  -N'p)    T'  (61) 

The  criteria  for  stability  specifies  that  all  the  C  . 
values  must  be  positive  and  that  the  following  determinant  is 
satisfied: 


I 
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0 


(02) 


0 


This  reduces  to  (C^  C^  -  C^  C^)  >  0  or 

Hence,  the  motion  will  be  stable  if 

-m/  (^^   +  d^)     n'   ^1  cTj  +  x;  r']>  rx/  y'  '^  -  /V  '^)  7''       (83) 

For  dynamically  stable  ships,  the  left-hand  side  of  the  inequal- 
ity is  positive;  consequently,  the  stability  or  instability  of 
the  towed  ships  is  essentially  a  matter  of  towpoint  location. 

For  dynamically  unstable  ships,  however,  the  extent  of 
instability  will  be  an  important  factor  in  determining  the  di- 
rectional stability,  since  the  sign  of  the  left-hand  side  of 
the  inequality  will  depend  on  the  sign  of  (i^  ^ \   ^2   '*^  ^ p     ^  ^■ 
Thus,  for  the  resulting  motion  to  be  stable,  (n^^  o^   +  x^   T  )>0, 
Hence,  the  dynamic  stability  criterion  ^1  must  be  small  if  this 
condition  is  to  be  satisfied. 

Reducing  towline  length,  however,  is  of  very  little  prac- 
tical consequence,  since  an  excessive  decrease  in  towline  length 
causes  hydrodynamic  interaction  between  the  towing  and  towed 
ships,  a  fact  that  will  cause  variation  in  the  hydrodynamic 
"constants".   Thus,  the  basic  assumptions  of  the  analysis  break 
down  and  the  conclusions  derived  become  untrustworthy. 
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NUMERICAL  ANALYSIS 

The  results  of  an  investigation  to  determine  the  hydrody- 
namic  coefficients  of  a  typical  barge  form  with  various  skeg 
configurations  made  on  the  rotating  arm  of  the  Experimental 
Towing  Tank  afford  the  data  necessary  for  an  analysis  based  on 
the  foregoing  treatment. 

The  particular  model  considered  had  been  found  to  be 
stable  according  to  the  free-towing  standards  set  by  the  Experi- 
mental Towing  Tank.  These  standards  stipulate  that  if  a  barge 
is  towed  from  its  bow  with  a  towline  length  of  twice  the  barge 
length,  that  the  barge  should  not  be  laterally  displaced  more 
than  -^  of  its  beam.  The  towing  device  is  a  carriage  that  pro- 
ceeds along  the  tank  rail  at  constant  speed  for  the  major  portion 
of  its  travel. 

The  accompanying  sketch  shows  the  bare  hull  of  the  barge 
model  together  with  the  skeg  modifications  investigated  on  the 
rotating  arm  (Ref.  6).   The  accompanying  tabulation  gives  de- 
tails of  the  configurations  as  well  as  the  comparable  results 
of  the  standard  free-towing  stability  criterion  set  by  the  Ex- 
perimental Towing  Tank  where  these  results  are  applicable. 


"T 


Cc 


II 


Case  I 


Case  111 


Case  IV 
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Model  Particulars: 
n'  =  .00219 
I  =  6.5  ft. 
P   =  1.9a 


«2  =  0.0351 

Xp  =  0.508 


■l  -  0.0219 

?'=  2.0 


Case  I   -  bare  hull  (unsatisfactory  by  free-towing  method) 

Case  II  -  bare  hull  and  skegs  (5°  to  ^  )  (unsatisfactory  by 
free-towing  method) 

Case  III  -  bare  hull,  skegs  and  stabilizing  prisms  (satisfac- 
tory by  free-towing  method) 

Case  IV  -  bare  hull  and  slotted  skegs  (not  tested  by  free- 
towing  method) 

Case  V   -  bare  hull  and  slotted  skegs  deflected  to  6  (not 
tested  by  free-towing  method) 

Case  VI  -  bare  hull  and  slotted  skegs  deflected  to  12°  (not 
tested  by  free-towing  method) 

In  each  case  the  model  was  tested  on  the  rotating  arm  at 

a  speed  length  ratio  of  ™=  =  0,1^63  where  V  is  measured  in 

knots.  The  following  values  of  the  hydrodynamic  coefficients 

were  determined: 


Case 

^''.^ 

N', 

0'=  t' 

>■'. 

V'; 

I 

0.006 

0.01;^' 

-0.00250 

0.0038 

-O.OO^:. 

II 

0.020 

0.0037 

-0.00280 

0.00U9 

-0.0029 

III 

0.020 

0.00U9 

-O.OO32U 

0.0055 

-O.OO3U 

IV 

0.022 

0.0069 

-0.00280 

0  0038 

-0.0028 

V 

0.020 

0.0060 

-0.00280 

0.0037 

-0.0029 

VT 

n.oiQ 

0.005 8 

-0.00280 

0.002Q 

-n.no^O 

Calculation  in  accordance  with  the  modified  development 
(pg.21)  and  Ref .  1  leads  to  the  following  values  for  the  untowed 
stability  criteria  (^i)  together  with  the  towed  stability 
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criteria  ^j  as  well  as  C .  and  Aj  . 


CASE    I 

CASE    l_l 

CASE    1 1 1 

CASE    1^ 

CASE  X 

CASE   Yi 

^1 

+  1   047 

+  0   237 

+  0  0735 

+  0   363 

+  0  303 

♦0   311 

^2 

2    IR*^ 

-2    131 

-2    195 

-2    268 

-2    197 

-2  222 

Co  X 

10^ 

0  76R7 

0   7«^fl7 

0   "'''^■' 

0.76P7 

0   7607 

0   76P7 

C,   X 

10^ 

0  8772 

1    4559 

1    6314 

1    4646 

1    455^^ 

1    46^*1 

Cp    X 

10^ 

-1    06 

1    270 

0.630 

0  023Q 

0    144» 

0    124P 

C^   X 

lO'^ 

-0.0532 

0.325Q 

0  4227 

0   3112 

0  3082 

0   2^34 

C4X 

10^ 

-0  0742 

0  0624 

0  Ofl52 

0  0602 

0  05«8 

0   0546 

10" 

0. 1426 

0   7029 

-0.2042 

-0.0261 

-0  0422 

Aj  X 

10'^ 

-0.0958 

0.0704 

-0.1926 

-0.  1326 

-0    1302 

^1 

0.82 

r   0.576 

r-0.0311 

r  0.  1254 

r     0.0905 

'     0  0919 

^2 

0.044 

J0.4775i 

L±0.509i 

[1  0 .  4606  i 

[+  0.  46341 

L  to. -45071 

^3 

*  0   24P  i 

-0. 1Q35 

-0.2405 

-0    1735 

-0.  1800 

-0.  1744 

■T 

4 

-1    873 

-1    8155 

-1    8195 

-  1    Q825 

-1    8^5 

-1    Q200 

The  tabulation  shows  that  none  of  the  cases  considered  is 
dynamically  stable  since  in  each  case  one  of  the  stability 
roots  (5r^)  is  positive.   Of  all  the  skeg  configurations,  only 
Case  III  produced  a  directionally  stable  barge  when  towed 
(since  this  is  the  only  condition  where  all  of  the  real  parts 
of  the  values  are  negative).   Case  III  was  also  considered 
stable  as  a  result  of  the  free-towing  tests  previously  des- 
cribed.  However,  the  drag  of  the  barge  in  this  case  is  notice- 
ably higher  than  for  any  of  the  others. 

It  will  be  seen  that  none  of  the  other  skeg  configura- 
tions creates  such  high  resistance.  While  these  other  cases 
are  not  directionally  stable,  it  appears  that  the  instability 
may  well  be  reduced  without  creating  excessive  drag.   It  does 
not  seem  too  unreasonable  to  assume  that  further  experimental 
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investigation  will  yield  skeg  forms  that  will  produce  direction- 
ally  stable  barges  without  adverse  drag  increments. 

The  presences  of  complex  "^^  values  (paired  in  the  table) 
indicates  that  the  motion  is  characterized  by  oscillatory  modes. 
In  all  cases  except  Case  III,  the  oscillation  is  undamped  and 
tends  to  build  up  with  time.   This  fact  imposes  continuous  al- 
ternate yawing  on  the  towed  ship.   In  Case  III,  which  is  di- 
rectionally  stable,  the  oscillations  still  appear  but  are 
damped  and  tend  to  diminish  with  time. 
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CONCLUSIONS 

The  preceding  analysis  shows  that  while  the  stability  cri- 
teria and  equation,  for  the  case  when  two  bodies  are  considered, 
have  been  determined,  their  evaluation  for  a  given  set  of  physi- 
cal and  hydrodynamic  constants  is  extremely  laborious.   It  is 
difficult,  if  not  impossible,  to  gain  an  extensive  amount  of 
qualitative  information  from  the  coefficients  of  the  stability 
equation  and  the  stability  criteria  because  of  the  complicated 
way  in  which  the  various  terms  comprising  these  expressions  ap- 
pear. 

In  the  simplified  analysis  where  only  the  towed  body  is 
considered  the  following  generalizations  may  be  made, 

1.  If  a  towed  ship  is  to  be  directionally  stable,  the  tow- 
point  must  be  located  forward  of  its  center  of  gravity  greater 
than  the  distance  between  the  center  of  gravity  and  the  center 

of  pressure  of  the  static  lateral  force. 

2,  Depending  upon  its  degree  of  stability,  mass,  towpoint, 
resistance,  and  moment  of  inertia,  a  dynamically  stable  ship  may 
become  unstable  when  towed  with  a  long  enough  towline. 

3o  Depending  upon  its  towpoint,  mass,  moment  of  inertia, 
and  damping  functions,  a  dynamically  unstable  ship  may  become 
directionally  stable  when  towed  with  a  short  enough  towline,  pro- 
viding that  there  is  negligible  hydrodynamic  interaction  between 
it  and  the  towing  ship. 

Ii.   The  motion  of  a  towed  body,  which  may  be  either  di- 
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rectionally  stable  or  unstable,  is  often  characterized  by  oscil- 
latory modes. 

5.  In  cases  where  a  ship  is  extremely  unstable  dynami- 
cally, no  possible  towing  configuration  can  make  it  direction- 
ally  stable . 

6.  Increase  in  hull  drag  aft  causes  variation  in  each  of 
the  hydrodynamic  parameters.   Sufficient  drag  will  cause  most 
towed  ships  to  become  direct ionally  stable.  The  experimental 
data  presented  indicate  that  the  drag  increments  may  be  mini- 
mized and  directional  stability  achieved  by  alteration  and  im- 
provement of  skeg  configurations. 
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APPENDIX  T 
NOMENCLATURE 

All  linear  dimensions  are  in  feet  and  angular  dimensions 
in  radians.  The  time  unit  is  the  second. 

J,       barge  length 

X    barge  towline  point  of  attachment  forward  of  barge 
center  of  gravity 

CG   barge  center  of  gravity 

Jf,        towboat  length 

Xp   towboat  towline  point  of  attachment  aft  of  towboat 
center  of  gravity 

CG  towboat  center  of  gravity 

0  barge  heading  angle 

^  barge  yaw  angle 

0  towboat  heading  angle 

/S  towboat  yaw  angle 

8    towboat  rudder  angle 
—  • 

0  towboat  rudder  angle  order 

t  time 

T  towboat  rudder  time  lag 

y  barge  linear  velocity 

V  towboat  linear  velocity 

1  towline  length 

L        over-all  length  of  tow,  i,e.,  from  bow  of  towboat  to 
stern  of  barge 

W  mass  of  barge 

W  mass  of  towboat 

Wj  longitudinal  virtual  mass  of  barge 

U  lateral  virtual  mass  of  barge 

W,  longitudinal  virtual  mass  of  towboat 

Wj  lateral  virtual  mass  of  towboat 

/  moment  of  inertia  of  barge 

0 

Tp   moment  of  inertia  of  towboat 

i    virtual  moment  of  inertia  of  barge 
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virtual  moment  of  inertia  of  towboat 

moment  of  inertia  of  displaced  fluid  of  barge 

moment  of  inertia  of  displaced  fluid  of  towboat 

longitudinal  virtual  mass  coefficient  of  barge 

lateral  virtual  mass  coefficient  of  barge 

longitudinal  virtual  mass  coefficient  of  towboat 

lateral  virtual  mass  coefficient  of  towboat 

virtual  inertia  coefficient  of  barge 

virtual  inertia  coefficient  of  towboat 

towline  tension 

density  of  water 

hydrodynamic  lateral  force  of  barge 

hydrodynamic  lateral  moment  of  barge 

hydrodynamic  lateral  force  of  towboat 

hydrodynamic  lateral  moment  of  towboat 

"■  first,  second,  and  third  derivatives  with  respect 
to  time 

^    indication  of  dimensionless  coefficient 

subscript  notation  denotes  a  partial  derivative 
with  respect  to  the  subscript;  thus 

dimensionless  longitudinal  virtual  mass  coefficient 
of  barge 

dimensionless  lateral  virtual  mass  coefficient  of 
barge 

dimensionless  longitudinal  virtual  mass  coefficient 
of  towboat 

dimensionless  lateral  virtual  mass  coefficient  of 
towboat 

dimensionless  virtual  moment  of  inertia  coefficient 
of  barge 

dimensionless  virtual  moment  of  inertia  coefficient 
of  towboat 


All  other  auxiliary  symbols  are  defined  where  they  appear 
in  the  text.  The  exact  expressions  for  the  dimensionless  factors 
are  given  in  each  instance  where  they  are  first  introduced. 
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